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' This paper studies the Schrodinger operator with Morse potential Vk{u) = \e^'^ + ke^ 

on a right half-hne [no,oo), and determines the Weyl asymptotics of eigenvalues for con- 
p stant boundary conditions at the endpoint uq. In consequence it obtains information on 

^ \ the location of zeros of the Whittaker function Wk^^{x), for fixed real parameters k, x 



Abstract 



with X > 0, viewed as an entire function of the complex variable /i. In this case all 
zeros lie on the imaginary axis, with the possible exception, if k > 0, of a finite num- 
ber of real zeros which lie in the interval — k < /i < k. We obtain an asymptotic for- 
mula for the number of zeros N{T) = {p \ Wk,p(x) = 0, \Im{p)\ < T} of the form 
N{T) = friogT + |(21og2 - 1 - logx)T + 0(1). Parallels are observed with zeros of 
the Riemann zeta function. 
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' 1- Introduction 

(N 

The one-dimensional potential function VA,Biu) = Ae^"^ + Be^ with ^ > and B real 
\ was proposed in 1929 by Philip Morse [SD] as an approximation to the (radial) quantum- 

mechanical potential for diatomic molecules, cf. Flugge |301 Problem 70]. From a quantum 
^ , physics viewpoint this corresponds to the Schrodinger operator 

viewed either on the line (— oo, oo) or on the left half-line (— oo, uq], with a suitable boundary 
condition imposed at the endpoint uq. By rescaling the independent variable we may 
without loss of generality reduce to the case A = j, and in that case relabel the other 
parameter, setting B = k, obtaining the potential Vk{u) = \e^^ + ke^. The Morse potential 
on the line, with A; = 0, arises elsewhere in physics in connection with the Liouville model 
on the hyperbolic plane (Grosche [34j) and in describing motion on the hyperbolic plane 
with a magnetic field (Ikeda and Matsumoto [35j). The latter authors also relate it to the 
Selberg trace formula on compact hyperbolic surfaces. 
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In both the hne and left half-hne cases the Schrodinger operator has an absolutely 
continuous spectrum of multiplicity one supported on the half-line [0,cxd), plus a finite 
(possibly empty) discrete spectrum, lying strictly below the absolutely continuous spectrum. 
On the line, the operator is essentially self-adjoint for the standard domain of smooth 
rapidly decreasing functions, and is exactly solvable. Formulas for its resolvent (Green's 
function) have been obtained, e.g. |35i Prop. 4.1(i)]. The absolutely continuous spectrum is 
described using generalized eigenfunctions given in terms of Whittaker Mfj,^(x)-functions, 
and the discrete spectrum has eigenfunctions given using rescaled Laguerre polynomials, cf. 
Ismail and Koelink [36] . 

1.1. Morse potential on right half-line 

We treat the Morse potential Vfc(n) on the right half-line [uo,oo), with constant boundary 
conditions imposed at the endpoint uq. The spectrum of this operator is quite different from 
the cases above: it is pure discrete, simple, and bounded below. This follows because the 
Morse potential is bounded below and is unbounded as u ^ cxd. The discrete eigenfunctions 
on the half- line are expressible in terms of Whittaker functions W^^^iix), in which /x is the 
eigenvalue parameter. We consider 



and in §2 we observe for all complex E there is a unique eigenfunction (up to scaling) in 
-L^([no, oo), dti) with eigenvalue E = — //^ given by 



The function 'ip{u,E) is well-defined since Wk^^{x) = Wk-fi{x) for all G C, when k is 
real and x = is positive real. The allowed real eigenvalues Eq < Ei < E2 < ■■■ in the 
discrete spectrum are selected by the boundary conditions at the left endpoint uq, with 
corresponding value xq = e"", since we are in the limit point case at the singular endpoint 
u = +00. 

This paper determines the Weyl asymptotics of the spectrum for fixed constant boundary 
conditions, obtaining in §4 that 



as T ^ 00. Dirichlet boundary conditions correspond to zeros of the Whittaker-function 
in the /i-variable, holding the other variables fixed. The results above, applied with these 
boundary conditions, lead to determination of the location and asymptotic distribution of 
the zeros of these Whittaker functions, given in §5 below. Namely, these zeros all lie on 
the imaginary axis, with finitely many exceptions, with asymptotics given in (jl.3p below. 
The exceptional zeros exist only if k = —k > 0, and they all lie on the real axis in the open 
interval (— |A;|). 




(1.1) 



(1.2) 



#{En <T} = ci\/f logT + csVf + 0(1) 
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1.2. Analogy with distribution of Riemann zeta zeros 

Our main interest in this family of operators arose from number theory, to view them as 
"toy models" for operators whose eigenvalues represent zeros of the Riemann zeta function. 
We note there has been a persistent effort to find a spectral interpretation of the Riemann 
zeta zeros. In this regard, by a "Hilbert-Polya" operator we will mean an (unbounded) 
self-adjoint operator on a Hilbert space whose spectrum encodes the zeta zeros, such that 
the self-adjointness of the operator encodes the Riemann hypothesis. Here we consider 
one-dimensional Schrodinger operators on a half-line from this viewpoint. 

To compare the distribution of Morse half-line eigenvalues with zeros of the Riemann 
zeta function C{s), in our context the s-variable in (^(s) corresponds to the /i-variable in the 
Whittaker function W,i^f^{x) via s := fi — ^; this variable change maps the imaginary axis 
in the //-variable to the critical line Re{s) = ^ in the s-variable. One point of this paper is 
that the function 

Zi{s) := Zi{s;k,xo) = W^,,^i{xo), 

for fixed real n and fixed real xq > has a number of properties analogous to the Riemann 
^-function 

^(,) i)^-fr(£)c(.). 

Namely, in Theorem 14. II we establish the following properties of Zi{s). 

(1) Zi(s) is an entire function of order 1 and maximal type, which is real on the real 
axis and is real on the line Re{s) = |. 

(2) Zi{s) satisfies the functional equation 

= Zi(l-s). 

(3) The number {Zi{s)) of zeros p of Zi{s) with \p\ < T is given as T ^ oo by 

Ar±(Zi(s)) = -TlogT + -(21og2- 1-Iogxo)r + O(l). (1.3) 
vr TT 

(4) All but finitely many of the zeros of Zi{s) lie on the line Re{s) = ^. The finite set 
of "exceptional zeros" all lie on the real axis, in the open interval (^ — ^ -|- |k|). There 
are no exceptional zeros if k < 0. 

(5) All zeros of Zi{s) are simple zeros, except for a possible double zero at s = ^. 

The zeros p = ^ + ij of Zi{s) = Zi{s; k, xq) have the following spectral interpretation: p is 
a zero of Zi{s; k, xq) if and only if = 7^ is an eigenvalue of the Schrodinger operator for 
the Morse potential Vfc(n) on the half-line [uq, 00) for k = —n and xq = e"°, with Dirichlet 
boundary conditions taken at the left endpoint uq. The self-adjointness of this problem leads 
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to the "Modified Riemann Hypothesis" property (4). This spectrum is bounded below by 
E > —K^, and the eigenvalues E < correspond to exceptional zeros of Zi{s; k,xo) lying 
on the real axis. 

In comparison with Zi{s), the Riemann ,^-function has properties (1) and (2). 
Analogous to property (3), the number of zeros N^{S^{s)) of the Riemann ,^-function with 
\p\ < T satisfies 

^t(C(^)) = -TlogT + i(-log27r-l)T + 0(logr), (1.4) 

vr TT 

cf. Edwards |261 Sec. 6.7] or Titchmarsh [68| Theorem 9.4]. Analogues of properties (4) 
and (5) are conjectured to hold for ^(s). Property (4) is a "modified Riemann Hypothesis", 
which permits a finite number of exceptional zeros to occur on the real axis. The Riemann 
.^-function has no exceptional zeros, but the possibility of exceptional zeros remains open 
for various generalizations of the Riemann zeta function, e. g. automorphic L-functions, 
cf. Iwaniec and Sarnak [37]. Concerning property (5), the zeta function is reported to have 
simple zeros for the first 10^^ zeros on the critical line. 

A possible weakness of the Morse potential as a "toy model" for a spectral interpretation 
of the zeta zeros is that its eigenvalue distribution only reproduces the main terms in the 
asymptotic number of zeros, it does not correctly model the "lower order" asymptotic 
behavior of zeta zeros. The detailed structure of scaled zeta zeros is conjectured to be 
described by the GUE distribution of random matrix theory, and the functions here do not 
reproduce this behavior. This detailed behavior presumably encodes all the number theory 
having to do with the distribution of prime numbers. This is discussed further in §5. A 
putative "Hilbert-Polya" operator Schrodinger potential, if it exists, will likely have exotic 
features to produce such features in its spectrum. 

1.3. Schrodinger operators and de Branges spaces as spectral models for 
zeta zeros 

Recently we observed that, provided the Riemann hypothesis holds, there is a natural can- 
didate for a "Hilbert-Polya" operator within the framework of the de Branges theory of 
Hilbert spaces of entire functions ([H], [IS]). The de Branges theory provides model op- 
erators for a certain class of operators, which includes generalizations of operators from 
singular Sturm-Liouville problems of the type considered here. This theory gives a rep- 
resentation of the associated operator as a (self-adjoint) canonical system of differential 
equations (defined in §5), in which each zeta zero p = ^ + corresponds to a simple real 
eigenvalue 7 of the resulting canonical system, so that the operator has a simple spectrum 
that is unbounded above and below. The results in (([41], [42] ) show that if appropriate 
Riemann hypotheses hold then the associated canonical systems must exist. We call these 
Hilbert-Polya canonical systems. 

The connection to Schrodinger operators is that, in favorable circumstances, a canonical 
system may be formally transformed to a Schrodinger operator on a half-line. This transform 
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squares the eigenvalues, so that the resulting Schrodinger spectrum is bounded below. If 
this transformation could legitimately be done for the putative "Hilbert-Polya" de Branges 
spaces above, then it would result in a "Hilbert-Polya" operator that is a Schrodinger 
operator on a half-line. If so, it would be very interesting to determine the associated 
potential of this operator. Note that this "Hilbert-Polya" potential, if it exists at all, may 
be highly singular, involve fractal features, and be describable only as a distribution. We 
remark that additional symmetries of the problem indicate that this Schrodinger operator 
necessarily would impose Dirichlet boundary conditions at the left endpoint. 

The Schrodinger operators treated in this paper may provide hints about features of this 
possible "Hilbert-Polya" potential. In particular, they give the dependence of eigenvalue 
asymptotics on the endpoint of the interval. Furthermore, even if such a "Hilbert-Polya" 
Schrodinger operator does not exist, these examples may still give hints about the canonical 
system of the de Branges space associated to RH. Namely, there exists a de Branges space 
with explicitly known canonical system for which the nonlinear transform above can be 
carried out to give a "toy model" Morse potential Schrodinger operator with k = ([43j). 
The form of this particular canonical system then gives hints about features of the associated 
canonical system coefficients of the Hilbert-Polya canonical system, which is guaranteed to 
exist if RH holds. 

1.4. Previous work 

There has been a long history of work on special functions whose zeros mimic the zeros 
of the Riemann zeta function, some of them associated to differential equations. Polya 
[55] observed in 1926 that the zeros of the iC-Bessel function K^{x) for fixed x > lie 
on the imaginary axis and have asymptotics similar to that of the Riemann zeta function. 
In fact Polya actually studies a function given by an integral representation, which in 
retrospect was shown to be a K-Bessel function, cf. Gasper [31]. In 1927 Polya [54] gave 
a second proof of the property that these zeros lie on the imaginary axis which made 
use of a second-order differential equation in an auxiliary variable. Another "toy model" 
discussed by Polya is described in Titchmarsh [68', Sect. 10.1]. It involves a shifted sum of 
two ii'-Bessel functions, discussed more recently in Gasper ^31j. The present work can be 
viewed as obtaining an extension to a one-parameter family of Polya's results, in the sense 
that a ii'-Bessel function is essentially a special case of the Whittaker functions considered 
here, namely K^{w) = ^y^Wo^^{2'w). Recent work of Biane [5] notes an operator-theoretic 
parallel between the i^T-Bessel function and the Riemann ^-function, and makes a very 
interesting new connection between these two functions in terms of probability theory. 

Aside from the "toy model" aspect, our results provide some new information on the 
location of zeros of Whittaker functions. Here there has been considerable previous work 
locating zeros of Wk_^(2;), in connection with various applied problems. This work mostly 
considers the case of fixed k and fi, and varies the (complex) variable z. In 1915 Milne [l9] 
studied real k and ^ and determined the nature of zeros in the z-variable on the positive 
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real axis in special cases. In 1941 Tvestkoff \71\ announced results on the location of all 
zeros with n and /i real and z = x > a positive real, and in 1950 Tricomi [69] gave 
detailed proofs. Their work determines the locations of the finite set of "exceptional zeros" 
in (4) above. In 1960 Dikii [21j and Dyson [25j independently considered the case of k real, 
fj, purely imaginary, and z a (nonzero) complex variable. Here Wi^^f^{z) is multi-valued in 
the z-variable with a singularity at z = 0, and the z-plane is cut along the nonpositive 
real axis. Their work arose from analysis of the stability of a hydrodynamic model of an 
incompressible atmosphere having density decreasing exponentially with height, with wind 
velocity increasing linearly with height, see Dikii [22] and Case [16] . repectively. Both Dikii 
and Dyson showed that all the zeros in the complex variable z in the region — vr < arg{z) < it 
lie on the positive real axis, and that there are infinitely many such zeros, having the singular 
point z = as a limit point. The results here complement these analyses by determining 
the location and density of zeros when k, is real and z is positive real. 

Our interest in de Branges spaces was influenced by recent work of Burnol (|10j. jllj . 
|13|). who described a scattering associated to the Fourier and Hankel transforms, and 
pointed out a de Branges space interpretation. Burnol constructed certain Dirac opera- 
tors and Schrodinger operators (corresponding to particular canonical systems) which have 
eigenvalue asymptotics of the general form cT log T + 0{T). In these papers Burnol did not 
explicitly note these asymptotics, but he elsewhere derived similar eigenvalue asymptotics 
for related operators ([12]). 

1.5. Methods 

The proofs given here are entirely classical, within the framework of Titchmarsh [67j . They 
break no new ground methodologically. The main point of this paper is rather to explain 
an analogy of these operators with a "Hilbert-Polya" operator for the Riemann zeta zeros. 
In particular, we explicitly determine the eigenvalue asymptotics (|1.3p . with all details in- 
cluded, for the convenience of number theorists. 

1.6. Contents of Paper 

The detailed contents of this paper are as follows. 

In §2 we determine the L^-eigenfunctions of the Schrodinger operator with the Morse 
potential on the half- line [tio,oo), and determine its solutions for constant boundary con- 
ditions at the left endpoint uq, particularly Dirichlet boundary conditions (Theorem 12. ip . 
This is a singular boundary value problem. Using Sturm-Liouville theory we deduce certain 
monotonicity properties of the eigenvalues as functions of the parameters (Theorem I2.2p . 

In §3 we show that the density of eigenvalues < T} of the Schrodinger operator is 
asymptotic to ciVTlogT + C2VT + 0{1), for certain constants ci,C2 (Theorem 13.21) . 

In §4, we specialize to Dirichlet boundary conditions, and apply these results to Whit- 
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taker functions Wk^^{x), viewed as a function of the complex variable fi. We establish that 
the zeros of Z{fi) := lie on the imaginary ^-axis, with finitely many exceptions, 

which themselves lie on the real axis, and we determine their asymptotic density to height 
T to be given by (|1.3p above (Theorem I4.ip . 

In §5 we discuss analogies of these results with zeta functions in number theory and 
with de Branges spaces. In this connection, it is well known that one can relate de Branges 
spaces to Schrodinger operators (Sturm-Liouville operators) in various ways, see Dym [23], 
Dym and McKean [23], and Remling [57]. Specific de Branges spaces having number theory 
content were related to Schrodinger and Dirac operators by Burnol [9], [10], [11], [13], [H]. 

In §6 we make some concluding remarks, discussing analytic spectral invariants associ- 
ated to these operators, the Weyl-Titchmarsh m-function of a Schrodinger operator and an 
analogous quantity in the de Branges theory, which we call the de Branges m-function. 

There are two appendices. Appendix A reviews properties of Whittaker functions 
Mk,^(z) and VFk,^(z) as functions of three complex variables {k, fj,, z). Appendix B gives for- 
mulas for the principal Weyl-Titchmarsh m- functions associated to Morse potentials Vk{u) 
on the right half line. This function encodes spectral data about the potential, and we 
show it is a ratio of Whittaker VF-functions with shifted K-parameter values, shifted by 
a constant. It is well known that a Schrodinger equation potential V{u) can be uniquely 
reconstructed from the principal Weyl-Titchmarsh m-function under general conditions. 

Notation. We follow the mathematical convention that Hilbert spaces have scalar prod- 
ucts < f,g > that are linear in the first factor and conjugate- linear in the second factor, 
as in Coddington and Levinson [T7] and Levitan and Sargsan [l5]. This is done for com- 
patibility with the theory of de Branges [3 p. 50] discussed in §5. References to Reed and 
Simon [HS]) who use the physics convention that Hilbert spaces are conjugate- linear in the 
first factor, must be adjusted for this fact. We set L'^{uq, oo) := L'^[{uQ,oo),du]. The vari- 
ables a, b, c, d, k, t, u, v, x, y denote real variables. The variables z = x + iy, s = a + it, k = 
a + ib, n = c + id denote complex variables. 

Acknowledgments. Christian Remling supplied helpful comments and references, and 
J.-F. Burnol made useful comments on his work and on m- functions. I thank them and the 
reviewers for helpful comments and corrections. 

2. Morse Potential on the Right Half Line 

Let k be an arbitrary real number. We consider the Schrodinger operator H = + 
for the Morse potential 

Vk{u) = ^e^- + ke'^ 

on the half-line [uq,oo). This is equivalent to treating the more general potential 

Va,b{^) = ^e^" + Be'' 
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with ^ > and B real on the half-hne [vq,oo), since by a translation of the time variable 
V = u — |log^, we reduce to a potential of the form Vk{u) above with k = on the 
half- line [uq, oo) with uq = vq — ^ log ^. 

Theorem 2.1. The Schrddinger equation 

-^ + Vk{u)H{u)=E^{u) (2.1) 

for the Morse potential Vfc(n) = ^e'^^+ke^ on a half-line [uq, oo) has the following properties. 

(1) The two-dimensional complex vector space £E(yk{u)) of solutions f{u) is given by 

(^(n)=e-t/(e") (2.2) 
in which f{x) is any solution to Whittaker's differential equation 

with parameters {k, jj,) given by k = —k and fi^ = —E. 

(2) For arbitrary E gC there is a one- dimensional subspace o/ iS'£;(Vfc(ti)) of solutions 
that belong to L^(no,oo). Setting E = , this subspace is spanned by a solution il){u,E) 
given by 

^P{u,z^) := e-tVF_fc,,,(e") = e~^W.k,-iz{e''), (2.4) 

where W^^^lx) denotes the Whittaker W -function. 

(3) The potential Vk{u) is of limit point type at +oo. If the constant boundary conditions 

(cos a)^/'(tio) + (sina)^/^ (uq) = 0. (2-5) 

for fixed < a < 27r are imposed at the left endpoint u = Uq, then the Schrddinger equation 
is self-adjoint. This operator has a pure discrete simple real spectrum, bounded below. For 
Dirichlet boundary conditions for parameter values k > 0, the smallest eigenvalue Eq > 0; 
for parameter values k < 0, one has Eq > —k'^. 

Remark. The case a = (resp. a = ^) gives Dirichlet (resp. Neumann) boundary 
conditions at the left endpoint. Recall that for general boundary conditions (j2.5p . there is 
no uniform lower bound on the spectrum as a is varied over < a < 27r. The eigenvalue 
Eq = Eo{a) is a continuous function of a except for a jump discontinuity at a = 0. 

Proof. (1) Whittaker's differential equation (j2.3p with parameters (k,/^) was studied in 
1904 by Whittaker [73]. Properties of its solutions are given in Whittaker and Watson [74^ 
Chap. 16], Buchholz [8]. Note that the parameters it// give the same solution space, since 
H only enters the differential equation via For fixed {K,fi) ()2.3p has a two-dimensional 
complex vector space of solutions SEiVkiu)) (for z on the positive real axis). This vector 
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space is spanned by the Whittaker functions (VFk^^(z), Mk^^(z)), except when k is a non- 
positive integer. (Another basis is {Wk~^{z),Mi^^-^{z)). Note that W^^^iz) = W^-fMiz) 
but in general Mf^^^^z) ^ Mk,_^(z). One calculates by direct substitution that if f{z) is 
any solution to Whittaker's equation ()2.3p . then ip{u) := e~2 /(e") satisfies the Schrodinger 
equation 

+ V-.iu)J <Piu) = -f^^cPiu). (2.6) 

Choosing (k, fi) = {—k, ±\/—E), these functions span the two-dimensional space of solutions 
SE{y-K{u)) to the Schrodinger equation (12. ip . 

(2) The Whittaker function Wf,^^{z) ([H Chap. 16]) is a solution f{z) of Whittaker's 
equation ()2.3p which has the property of rapid decay as z = x — > oo (on the real axis), 
made unique by the normalization of its asymptotics 

W^«,^(x) = e-5-x'^(l + 0(^i)) (2.7) 

as X — > oo. Viewed in the complex domain, the same asymptotics holds for \arg(z)\ < ir—e < 

TT, with the implied constant in the error term depending on e > 0. All other linearly 

independent solutions to Whittaker's differential equation increase exponentially as x ^ 

oo. Choosing k = —k, the particular choice ip{u) := 'ip(u,z'^) coming from the Whittaker 

function is the only solution (up to a scalar multiple) that satisfies ^{u) G L^(no,oo) for 

1 ^ 

(any) finite uq; it has double-exponential decay as n — > oo, roughly like e~2'^ , as specified 
by the asymptotic formula (j2.7p . All other solutions have double-exponential growth as 
u — > oo. The solutions (j2.4p thus parametrize the subordinate solutions of the Schrodinger 
operator with Morse potential on the right half-line, in the sense of Gilbert and Pearson 



(3) The property that Vfc(n) is positive for sufficiently large u puts the Schrodinger 
equation in the limit point case at n = oo (cf. Reed and Simon [55, Theorem X.8]). In 
particular for each real E there is at most one solution to the Schrodinger equation with 
eigenvalue E (up to a scalar multiple) that belongs to L^([uo, oo); du). Here such a solution 
was exhibited for all complex E, given above. 

It is known that for continuous potentials V{u) on [uq, oo) that are bounded below, the 
Schrodinger operator + V{u) with (separated) boundary conditions given by 

(cos a)f{uo) + (sin a)/ (uo) = (2.8) 

for fixed < a < 27r is self-adjoint in L'^{[uq, oo); du). To describe the explicit domain of this 
operator, let ^C[uo,oo] denote the set of absolutely continuous functions f{u) on [uo,oo) 
whose (almost everywhere defined) derivative is in L'^[uo, oo). The associated dense domain 
for self-adjointness is 

Va := {f{u) £C'^[uo,oo)nL'^{[uo,oo);du) \ f'{u) £ AC[uo,oo], 



+ V{u) ] f{u) G L^{[uo, oo); dn), and i^M) holds}. 
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cf. Reed and Simon [551 P- 144]. 

The fact that Vk{u) is continuous, with Vk{u) — > oo as u — > oo is sufficient to imply 
that it has pure discrete spectrum, bounded below, for the boundary conditions ()2.5p . for 
any choice of a (Levitan and Sargsan |45[ Sect. 4.1, Lemma 1.2].) This spectrum is simple 
because for each real E there is at most one solution in L'^{[uo,oo);du). The spectrum is 
bounded below, for any fixed constant boundary conditions, because the potential Vk{u) is 
bounded below. 

The Sturm-Liouville theory applies to real potentials having a pure discrete simple 
spectrum bounded below. It implies that the eigenvalues can be numbered Eq < Ei < 
E2 < ... and the corresponding eigenfunctions ipo < ^pl < ■02 < ••• can be chosen real- 
valued. The Sturm-Liouville theory shows that the eigenfunction ipn{u) has exactly n sign 
changes on the domain [no,oo), so has n zeros; the bottom eigenfunction ipoiu) has no 
changes of sign and can be chosen positive. (See Atkinson ^ Chap. 8]), and Coddington 
and Levinson [17, Chap. 8].) 

We now impose Dirichlet boundary conditions q = 0, and lower bound the minimum 
eigenvalue Eq using the Sturm comparison theorem. For A; > the eigenvalues of Vk{u) 
are bounded below by those of Vo{u) (with same boundary conditions on the endpoint uq), 
since Vk{u) > Vq{u) on the entire line. It suffices to show that the minimum Dirichlet 
eigenvalue E > —k"^ for Vfc(n) with k < 0, and Eq > when A; = 0. The function Vk{u) 
for k < has global minimum value —k"^, which occurs at the unique point u = log(— 2fe), 
while for k = 0, Vk{u) is strictly positive. We use the Sturm comparison theorem for the 
potential Va,uo,e{u) '■= ^{u — ^0)^ — S for ^ > 0, so that — ^ + Va,uo{'^) on [no,oo) is 
just a translated and scaled version of the harmonic oscillator potential on [0,oo), with 
eigenvalues shifted by E. For a given e > , we take E = k"^ + e and then can choose ^ > 
just small enough that 

Va,uo,e{u) = A{u - uof - k^ -e< Vk{u), uq < u < 00. 

The Sturm comparison theorem in the singular case (see Coddington and Levinson [17^ 
Chap. 8, Theorem 1]) now shows that the minimum eigenvalue Eq of Vfc(n) has 

EQ>EQ-k^- €, (2.9) 

where Eq is the minimum eigenvalue of VAfl,o{u) = Au^ on [0, 00) with the same boundary 
conditions at n = liQ. It is well known that the oscillator Hamiltonian + on the 

half-line [0, 00) with Dirichlet conditions at n = has a strictly positive spectrum with 
minimum eigenvalue Eq = 3. Under rescaling of A we retain the property Eq > 0, and the 
inequality Eq > —k'^ then follows from (j2.9p on letting e become sufficiently small. Finally, 
for A; = we may omit e, and so obtain Eq > 0. ■ 

The Sturm theory implies monotonicity properties of the eigenvalues as the parameters 
k, Uq are varied, with boundary conditions remaining fixed. 
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Theorem 2.2. Let Vk{u) = je^" + /ce". Let En^a{uo,k) denote the n-th eigenvalue (n > 0) 
from the bottom of the spectrum of the operator H = —-^ + Vfe(u) on the half-line [iio,oo) 
with boundary conditions 



with a being fixed. For each n > the following hold. 

(i) For fixed parameter k > 0, the n-th eigenvalue En^a{uo,k) is a strictly increasing 
function of uq . 

(a) For fixed parameter k < 0, the n-th eigenvalue En,aiuo,k) is a strictly increasing 
function of uq, when uq > \og2\k\. 

(Hi) For variable k, with fixed uq, the n-th eigenvalue En^a{uQ, k) is a strictly increasing 
function of k. 

Proof. By a change of parameter v = u — uq we shift the left hand endpoint to v = 0. 
Then each of the assertions (i)-(iii) above can be viewed as varying the potential, while 
keeping the left endpoint and the boundary condition fixed. In all cases the potentials are 
bounded below, and the variation increases the potential pointwise. The results then follow 
from comparison of the minimax characterization of the n-th eigenvalue. (Reed and Simon 
[Ml XIII. 1]). □ 

3. Asymptotic Formula for Eigenvalue Density 

We obtain estimates for the number of eigenvalues using a version of Weyl asymptotics, 
which says that the number of eigenvalues < T of a Hamiltonian H = p^ + V{q) with 
V{q) bounded below can be approximated by the volume in the phase space of the region 



The half-line condition requires that we restrict the g-variable in the phase space to 
[no, oo) in estimating the phase space volume. The following result is a slight modification 
of a result of Titchmarsh [67| Theorem 7.4]) giving a rigorous estimate of this kind. 

Theorem 3.1. Let V[u) be a real-valued C"^ -potential on [uq^co) which for all sufficiently 
large u is increasing and convex downwards. Then the Schrddinger operator —-^ + V{u) 
on [uo,cxd) with boundary condition 

(cos Q)V'(^to) + (sinQ)'i/''(uo) = 0, 

has pure discrete simple spectrum. Let N{T; a, uq) count the number of eigenvalues En < T 
for this boundary condition. Then there is a constant Ti such that for all T >Ti there is a 




(2.10) 



H <T. 
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unique solution ut > uq to V{u, T) = T, and 

1 

N{T;a,uo) = - ^/T - V{u)du + 0(1) (3.1) 

holds for Ti < T < oo. 

Proof. Suppose that V{u) is increasing and convex downwards for u > ui. The convexity 
condition imphes for u > ui that V{u) is strictly increasing and increases at least at a 
linear rate. Thus V{u) oo as u ^ oo, so the Schrodinger operator with given boundary 
conditions has a pure discrete simple spectrum. 

Titchmarsh \67\ Theorem 7.4] proves the estimate ()3.ip assuming as hypotheses that 
V{u) is continuous, increasing and convex downward on the whole interval [uq,oo), i.e. 
that one can take ui = uq. A slight modification of his argument handles the case above. 
We set Ti := max{y(n) : uq < u < ui} For T > Ti, the strict increasing property of V{u) 
for u > ui implies that the equation V{u) = T has a unique solution, so that the set of 
values {u : V{u) < T} is an interval [uojUt]- Titchmarsh's argument estimates the number 
of oscillations in the real- valued -solution iIj{u;E) with eigenvalue ii^ on [no, oo) (not nec- 
essarily satisfying the boundary conditions). We replace his integral estimate < /2(^) ^ | 
with hiY) = 0(1); the extra C^'^-condition on V{u) on the interval [uo,tti] is imposed to 
guarantee that the relevant integrals exist. ■ 

We apply this result to the Morse potential. Note that the parameter k in the Morse 
potential has no influence on the first two terms in the resulting asymptotic formula (j3.2p 
below. 

Theorem 3.2. For the Morse potential Vk{u) = |e^" -|- ke^ on [tio,oo) with boundary 
conditions {cos a)^p{uo) + {sin a)ip' (uq) = 0, the eigenvalue density satisfies 

N{T; a, uq) = -VTlog Vf + - {2log2 - 1 - uq) Vt + 0(1). (3.2) 
vr vr 

as T ^ oo, and the implied constant in the O- symbol depends on both k and uq. 

Proof. The Morse potential Vk{u) is increasing and convex downwards for all sufficiently 
large u, so Theorem 13.11 applies to give 

N{T, a) = - ^/T-Vk{u)du + 0(1). (3.3) 

^ Juo 

as T ^ OO. For k > the Morse potential is increasing and convex downward on the entire 
real line. However when A; < there is a critical value ui to the left of which this no longer 
holds. 

It remains to estimate the integral in (j3.3p . with Vk{u) = |e^" -|- fee". We first suppose 
that Uq is such that Vk{u) is increasing on [tio,oo), which holds for all uq if /c > and for 
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■"0 > 1*1 '■= log2|/c| if k < 0. We assume in what follows that uq > log(6|fc| + 1), and treat 
the case when uq < log(6| A:| + 1) at the end of the proof. We suppose T large enough to be 
in the region where Vfe(n) = T has a unique solution and introduce the positive constant 
T* given by (r*)^ = T + A;^ . We evaluate the integral 



I := ^T-Vk{u)du = £^ J (r*)2 - Qe« + 



2 

du. (3.4) 



'Uo JUQ 

We change variable u io w determined by 

w = 



^e^ + k 



Since Vk{u) is monotone increasing on [uq,oo) the change of variable is single- valued, and 



on the integration interval w decreases from wq = i ^ to wt = 1. We have 



and obtain 



(r*)2 / 1 \ dw 

w'^ 1 1- ^ j V 



Wo ' \ T' 



wo 



= T* X 1- 



1 1 \ dzi; 



1 V ■u^^yl-^/'u; 
in which — = ^e'^° + k. We let = e'', and obtain 

wo 2, ' 



/ = T*y^ 713^ j^^-^ j d^, (3.5) 
in which 

t-o = logw;o = logT* - log(^e"" + k), (3.6) 



where by hypothesis ie"° + > + i. Now 



2^ 



so 

1 



ie" + A; 



2=" + tje<'. 

The hypothesis on -uq > log(6|A;| + l)when A; < ensures that 



0<^<2- (3-7) 
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Thus we may legally expand the integrand of (|3.5p in a power series in ^ and obtain 
I = M + R with 



Vf := T* (^j^ ° Vl - e-'^^dv 



Now we have 



M = T*vo + T* (^j^ - - l)dv^ 

= T* (^logT* -log(ie"» + 

in which the constant in the O-symbol depends on k and uq. In the last line we used the 
estimate 

1 _o„ . . „ / 1 



/"OO /"OO -| 

/ {Vl^e^ - l)dv\ < -e-^''dv<0 

Jvn J\oRT*-ca 2 



/do ^logT*-co 

Now we use the identity 



*\2 



/ (1 - \/l - e-2^)(i^; = 1 - log 2, 

JO 



to obtain 

M = T*logT* + (^Iog2-1-Iog(^e"o + A;)^T*+O(l), (3.8) 

and the constant in the O-symbol depends on k and uq. The remaining term is split as 
i? = i?i + i?2, in which 

^1 - T'ifii^y r^'-dv 



The first term above is 



= T*(-log(l-^)+log(l-A)) 

= T* (-log(l-^-^ ) + 0{^)] 

= T* ('log(ie"'' +k)- log(^e"'')) + O (1) 
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The second term is estimated by 



i-2 



< o{i) + ^y( , y =0(1). 

Putting these estimates together gives 

R = R,+R^ = T* (log(ie"« + A;) - log(ie««)) + 0(1), 

and 

I = T* log r* + (2 log 2 - 1 - uo) T* + O (1) , (3.9) 

provided that uq > log(6|A;| + 1). 

In the remaining case where uq < log(6|A;| + 1), we split the integral into a portion from 
uq to ui = log(6|A;| + 1), and from [ni,oo). The integral from [ui,oo) is estimated by (|3.9p 
while that on [iio,Mi] by 

/ ' y/T-Vkiu)du = r yjT-0{l)du = {ui - Uo)Vf + O ( ^) , (3.10) 
Juo Juo WTJ 

in which the implied constant in the 0-symbol depends on uq and k. Note that as T* — > oo, 
we have 

T* = V¥T^=Vt(^1 + ^ + 0{^)^ =VT + o(^-^y (3.11) 

Combining this with (I3.10p shows that ()3.9p remains valid for all intervals [uqjOo). Next, 
substituting([33I]) in ^3M yields 

/= log\/r + (21og2-l-uo)\/r + 0(l). (3.12) 

Substituting this estimate for (j3.4p in (j3.3p gives the desired estimate (|3.2p of N{T, a). ■ 

4. Zeros of Whittaker Functions in the yU- Variable 

For Dirichlet boundary conditions the eigenvalues of the Morse potential on the half-line 
are specified by Theorem 12. II in terms of the zeros of Whittaker functions. We combine this 
with the results of §4 to deduce the distribution of zeros of these Whittaker functions in 
the ^-variable. The function Zi{s) in the introduction is given as Zi{s) = Z{s + ^), where 
Zifi) := VF«,^(e"«) below. 
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Theorem 4.1. For fixed real parameters k and uq the Whittaker function 

Z(/i) := W^,^{e^°) (4.1) 

is an entire function of /x of order 1 and maximal type. It is real-valued on both the real axis 
and imaginary axis. It is an even function of n, i.e. it satisfies the "functional equation" 

Z{^) = Z{-fi). (4.2) 

(1) The zeros of Z{jjL) lie on the real and imaginary axes. All zeros of Z{fi) are simple 
zeros, except for a possible zero at fi = 0, which if it occurs will be a double zero. 

(2) For each k the function Z[ji) has finitely many real zeros. There are no real zeros 
when K < 0, and for k >0 they all lie in the open interval —k < fi < k. 

(3) There are infinitely many imaginary zeros. The total number N{Z{fi)) of zeros with 
\Im{fj,)\ < T satisfies the asymptotic formula 

N{Z{fi)) = -riogr + -(21og2-l-no)T + 0(l), (4.3) 
vr vr 

as T ^ oo, where the O- constant depends on k. 



Proof. We set x = e^. It is well known that for fixed k and x the Whittaker func- 
tion Wfi,^{x) is an entire function of /j,, see Appendix A. The function Wk-^{x) satisfies 
Whittaker 's differential equation with the same parameters (k, /x), and it also has the same 
asymptotics as x — > oo, namely 



W^,^{x)=e-'^^x^{l + o(^^y, 



It follows that 

W^,^^{x) = W^,^{x) (4.4) 

holds for all parameters (k, ^) and positive real x. (It then holds universally under analytic 
continuation in the x-variable, see Appendix A.) The functional equation (j4.2|) follows, 
taking x = e^° . This also implies that the function G{fi) := Zi^^fji) is an entire function of 
/i. 

We next establish the real symmetry when k is real-valued and x = e"" is positive-real 
valued (on the principal branch of the Whittaker function). For 2[i not an integer we have 
the representation ([74l Sec. 16.41]) 

W.^^{x) = -^t^±-M,^^(x) + , M.,_^(x). (4.5) 



Now we assert that, for real k and positive real x, applying complex conjugation to ()4.5p 
yields 

W,,-^{x) = W^A^. (4.6) 
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To justify ()4.6p , we note that for positive real x, on the principal branch of the Whittaker 
M-function we have 

(see ()7.5p in the Appendix) and that T{s) = T{s) for s G C, since it is real on the real 
axis. This justifies (|4.6p . and in fact for real k one has W^^f^iz) = Wk,p{z) for complex z 
with I arg(2)| < vr. We conclude using (14. 5p that for real k and positive real x, the function 
Wk^^{x) is real when /i is on the real axis. In addition, under the same assumptions, the 
conjugation symmetry (j4.5p together with the symmetry ()4.4p implies that WK^^ix) is real 
when fi = it is on the imaginary axis. 

For real k and positive real x, a growth bound deduced from the contour integral repre- 
sentation given in formula ()7.9p of Appendix A implies that the function Z{n) = WK^^{e^°) 
is an entire function of order at most 1. The fact that Z{fi) is entire of order 1 and maximal 
type will follow from the zero counting formula (|4.3p . asserting that for some C > it has 
at least CTlogT zeros in the disk \z\ < T as T ^ oo. For any entire function of order 
less than 1 or of order 1 and finite type necessarily has 0(T) zeros in each disk \z\ < T as 
r ^ oo, (cf. Titchmarsh ^ 8.75]). We now establish (l)-(3). 

(1) Apply Theorem 12.11 for the Morse potential Vk{u) with k = —n. Then the zeros 
of Z{pi) := WK^^(e"") correspond to the L^-eigenfunction tp{uQ, —^J^) = 0, corresponding to 
E = — /i^ being a Dirichlet boundary condition eigenvalue of the Morse potential Schrodinger 
equation on the half-line [uq, oo). This boundary value problem is self-adjoint, so E is real. 
Thus n must be real or pure imaginary. 

The Schrodinger operator spectrum is simple, therefore each value of E occurs once. 
Both roots in ^ of -E = — produce a zero of the Whittaker function, due to the symmetry 
Wk,i_i{x) = Wf^-^{x). These zeros must therefore be simple, except for /u = 0, where a 
double zero must occur since Wk-^{x) is an even function of 

(2) Theorem 12.11 establishes that the Schrodinger spectrum for Dirichlet boundary con- 
ditions is pure discrete and is bounded below hy E = —fi^ > — A;^, and by -E > when 
A; > 0. Thus there can only be finitely many roots /i on the real axis, all with < fi"^ < k^, 
so that — |k| < < and there are no real roots when A; > 0, i.e k < 0. 

(3) The asymptotic estimate (|4.3p follows from Theorem 13.21 A pair of zeros fi = ±it 
corresponds to a single Dirichlet eigenvalue E = —fi'^ = t'^ oi the Schrodinger operator, 
and the condition < T corresponds to the bound E < T^. Thus we obtain from the 
eigenvalue asymptotic p.2p that 

iV(Z(/i)) = N{T; 0, no) = -Tlog T + -(2 log 2 - 1 - uo)T + 0(1), 

TT TT 

as required. ■ 
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Remark. For the case A; = 0, and ^(w) > the Whittaker function is related to the 
-RT-Bessel function (MacDonald function) by 

Kf.H = ^^oA"^^)- (4-7) 

The X-Bessel function for positive real w = is known to have all its zeros on the imaginary 
axis. According to Erdelyi et al, [281 Vol II, 7.13.2 (19)] an asymptotic formula is available 
when ^ = it IS on the imaginary axis with w = x = e^. This formula states that, for 
t > a; > 0, 

Ku{x) = V2^{t' - xY'/'e--2-' fsin(t cosh-i(-) - Vt^-x^ + ^) + 0( )) . 

\ X 4 y/t^ -X^ J 

(4.8) 

(We have changed the error term, which appears to have a misprint.) These asymptotics 
agree with that of Theorem 14.11 

5. Eigenvalue Interpretations for Riemann Zeta Zeros 

There is considerable circumstantial evidence for the existence of a spectral interpretation of 
the Riemann zeta zeros, see Berry and Keating [2] and Katz and Sarnak [38]. Indeed there 
are several known operator-theoretic interpretations of zeta zeros. There is a scattering 
theory interpretation of the Riemann zeta zeros arising from work of Pavlov and Faddeev 
|52| concerning the Laplacian acting on the modular surface. A slightly different version of 
this interpretation is given by Lax and Phillips [H], see in particular [HI Appendix 2 to Sect. 
7]. Recently an operator formulation of the Lax-Phillips approach was given by Y. Uetake 
[72] . There are also interpretations of the zeta zeros in terms of eigenvalues of operators 
on various function spaces, given by Connes [18]. The formulation of Connes permits the 
construction of such an operator on a Hilbert space that has eigenvalues that detects only 
those zeta zeros that are on the critical line, up to a fixed multiplicity. This approach 
was extended to automorphic L- functions by Soule [65] and Deitmar [20j . R. Meyer ([47], 
[48j ) gave an unconditional formulation of an operator on a more general Banach space 
whose eigenvalues detect all zeta zeros, including those that are off the critical line if the 
Riemann hypothesis fails. There is also recent work of Sierra ([SU]- [52]) and Sierra and 
Townsend [63] concerning quantum mechanical models for the zeta zeros, and possible 
"Hilbert-Polya" operators. Burnol [15] has constructed a specific "toy model" operator of 
the integral-differential form considered by Sierra, based on his earlier work [13j . 

We recently observed ([S], [13) that there is a natural candidate for a "Hilbert-Polya" 
operator, using the framework of the de Branges Hilbert spaces of entire functions, provided 
that the Riemann hypothesis holds. This interpretation leads to a possible connection with 
Schrodinger operators on a half line, described below. 

We first review some aspects of the de Branges theory. The de Branges theory of Hilbert 
spaces of entire functions ([7|) assigns to an entire function E{z) of a special type called a 
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structure function a reproducing kernel Hilbert space of entire functions 7i{E{z)), together 
with a (generahy unbounded) hnear operator : f{z) zf{z) defined on the domain 
Vz := {f{z) e Ti.{E{z)) : zf{z) G 7i{E{z))}. This operator is closed and symmetric, with 
deficiency indices (1,1), so that it has a one-parameter family of self-adjoint extensions, 
parametrized by the unit circle = U{1). A structure function E{z) is any entire function 
satisfying 

\E{z)\ > \E{z)\, when Im{z) > 0. (5.1) 

Functions E{z) with this property are sometimes called Hermite-Biehler functions (see Levin 
Chap. VII]). Any entire function can be uniquely represented as E[z) = A{z)—iB{z), in 
which A{z)^B{z) are entire functions that are real on the real axis. If we let E^{z) 
then 



Eiz) 



A{z) = \{e{z)+E\z) 



B{z) = ^.{e{z)-E\z) 



The Hermite-Biehler property (jS.ip has the consequence that the associated A{z) and B[z) 
have only real zeros, and these zeros interlace (allowing multiple zeros, counting multiplic- 
ity) (cf. [ini Lemma 2.2]). The Hilbert space TC{E{z)) is always nonempty, and it consists 
of all entire functions satisfying certain growth restrictions (compared to E{z), E'^{z)) in 
the upper half plane. One particular self-adjoint extension of (M^,!?^) has pure discrete 
simple spectrum located at the zeros of A{z), and another has pure discrete simple spectrum 
located at the zeros of B{z). 

The de Branges theory also gives a Fourier-like transform, depending on the structure 
function, that gives an isometry to a new Hilbert space on which the multiplication operator 
becomes a 2 x 2 matrix system of linear differential operators, called a canonical system; 
we call this the de Branges transform. Here a canonical system is a family of differential 
equations, depending on z G C as a parameter, on an interval [a, b] in which we allow 
< a < 6 < oo, given by 



in which 



di 


1 



Ait,z) 
B{t,z) 



-1 




z3M{t) 



M(t) 



A{t,z) 
B{t,z) 



a{t) (3{t) 
m lit) 



(5.2) 



with M(t) being a positive semidefinite real matrix- valued function on the interval {a,b), 
with coefficients being measurable real valued functions of i. In fact de Branges 's formu- 
lation of his theory [7] formulates results using an integral equation which is an integrated 
form of the canonical system, in order to deal with smoothness issues; he does not use canon- 
ical systems. However the use of a canonical system makes comparison with Schrodinger 
operator formulation easier, cf. Remling [57]. In general the de Branges transform to the 
canonical system is not known explicitly, but it has been determined in a number of "exactly 
solvable" examples, see de Branges [6], (Tj Chap. 3]. 
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A special subclass of canonical systems can be nonlinearly transformed to a pair of 
Schrodinger operators on an interval or a half-line, a transformation that squares the eigen- 
values. Such canonical systems are a subclass of those whose canonical matrix M(t) is 
diagonal and invertible almost everywhere, de Branges [71 Section 4] proves a result stating 
that a structure function has an associated canonical system that is diagonal whenever A(z) 
is an even function and B{z) is an odd function. For such diagonal canonical systems, one 
may monotonically rescale the time variable so that the coefficient matrix becomes 



M(t) 



ait) 
7(0 



(5.3) 



with determinant 1 almost everywhere, on a rescaled interval [d,b] with — oo < d <b < oo. 
If these coefficients are smooth enough, then the transform gives a (supersymmetric) pair 
of transformed Schrodinger operators — ^ + y^{t) having potentials 



with associated superpotential 

Wit) 

d 



1 a' it) 



2 ait) ' 

where a'(t) = ^a(t). These two Schrodinger operators have related boundary conditions at 
the left endpoint, such that they are self-adjoint with identical discrete spectrum, with the 
exception of the eigenvalue 0. One of them has Dirichlet boundary conditions, the other has 
boundary conditions that depend on the eigenvalue. There are a number of known examples 
where, starting from the structure function, all these steps can be carried out explicitly to 
obtain such a pair of Schrodinger operators. This is the case for the Morse potential for 
certain parameter values, using the fact that it is a shape-invariant potential treatable by 
supersymmetry, cf. Cooper et al. [121 Sec. 4.1], see [13]. 

In [10] , [H] , [12] we considered de Branges spaces associated to the Riemann zeta func- 
tion and more general L- functions. In particular, we associate to the Riemann ^-function a 
family of entire functions depending on the real parameter h, given for h ^ hy 

Ehiz) := \ [ii\ +h-iz)+ii\-h- iz)^ - ^ [ii\ + h-iz)-ii\-h- iz] 

and for /i = by 

E^{z) :=C(^-iz)+C'(^-iz), (5.5) 

with ^'(s) = -^^is). We show in ([40l Theorem 2.1]) that for > ^ each member of 
this family is unconditionally a de Branges structure function, and, conditionally on the 
Riemann hypothesis, that they are structure functions for all nonzero real h. Furthermore 
we show in ([HI Theorem 1]) that the Riemann hypothesis holds if and only if the entire 
function E^iz) in (j5.5p is a de Branges structure function. Assuming RH, a theorem of de 
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Branges (^Ti Section 47]) applies to assert that the canonical system associated to Eq{z) by 
the de Branges transform has matrix function M(t) of diagonal form. 

To summarize, if the Riemann hypothesis holds, the de Branges theory predicts the 
existence of a "Hilbert-Polya" operator given as a diagonal canonical system, specifically 
associated to the structure function (jS.Sp . Now we can ask if one can nonlinearly transform 
this diagonal canonical system to a pair of Schrodinger operators. As obstacles, we do 
not know whether the resulting diagonal matrix function M(t) will be almost everywhere 
invertible, and even if it is, we do not know whether the resulting matrix coefficients ()5.3p 
will be smooth enough to permit the nonlinear transform ()4.5p . Here we speculate that 
the smoothness conditions will be the main obstacle, and that one may wish to consider 
this transform in the sense of generalized functions. (This is a suitable topic for further 
research in its own right.) In any case, if this transform were possible, there would then 
exist a "natural" Schrodinger operator on a half-line, having a (possibly distributional) 
"Hilbert-Polya" potential that, for Dirichlet boundary conditions, has as simple eigenvalues 
7^ corresponding to p = ^ ±27 being a pair of zeta zeros. One can then attempt to reverse- 
engineer this potential, or in any case, to reverse-engineer the matrix coefficient functions 
of the associated canonical system. 

For the Morse potential treated here the "toy model" aspect can be extended further, 
to the level of de Branges spaces and canonical systems. That is, there is a family of de 
Branges spaces associated to the Morse potentials, described in [3^, where the nonlinear 
transformation to a Schrodinger operator on a half-line can be justified. As an example, 
one can show that for real u the function 

Eo{u, z) := e5('="-")vri_.^(e") - ize-^('^"+"%_ 1 ^,^(e") (5.6) 

is a de Branges structure function, and the transforms above convert it to a Schrodinger 
operator with Morse potential with parameter /c = — ^ on a the corresponding half-line 
[u, 00). In particular, the self-adjoint extension of the de Branges multiplication operator 
associated to 

A{z) = e-i('="-")Tyi_.^(e") 

is transformed to this Schrodinger operator with Dirichlet boundary conditions at the left 
endpoint. This example shows that there is no apparent obstacle at the level of asymptotic 
distribution of eigenvalues to finding a diagonal canonical system transformable to an as- 
sociated Schrodinger potential that encodes the zeta zeros. Such an (integrated) potential 
must increase exponentially as t ^ 00 to have the correct asymptotics of zero densities. 

We next address the feature that these "toy models" reproduce only the main term 
asymptotics of the distribution of zeta zeros, but not the fine structure of these zeros. 
Namely, the Riemann zeta zeros are believed to have local density statistics with fiuctuations 
described by the GUE distribution of random matrix theory, see Katz and Sarnak [38] . 
These GUE fluctuations encode properties of the distribution of primes, and they partially 
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manifest themselves in the fact that the error term O(logT) for zeta zero asymptotics in 
()1.4p is nontrivial. In contrast the zeros of the Whittaker functions have normalized local 
density statistics which provably have no fluctuations, i.e. the normalized zeros will be 
completely regularly spaced, a fact which corresponds to the 0(1) error term in (jl.3p . This 
fact is not so far removed from the zeta function itself, because in [41j a similar phenomenon 
was observed for zeros of differenced Dirichlet L-functions, whenever one moves off the 
critical line to a line Re{s) = ^ + 6, for some S > 0. The zeros of these differenced L- 
functions correspond to vanishing of the real part of the L- function + S + it,x)- The 
asymptotics of these differenced L-function zeros is known to have the same main term 
as the zeta zero asymptotics, but also to have a 0(1) remainder term, just as in these 
"toy models". Furthermore 01] showed that these differenced -L-function zeros do have 
an eigenvalue interpretation for a multiplication operator in a suitable de Branges space, 
assuming RH. This fact shows that the GUE phenomenon for the zeta zeros is (in some 
sense) confined to the critical line Re{s) = |. We take this to mean that, assuming RH 
holds, the de Branges space with structure function (jS.Sp will have some unusual features of 
its associated canonical system. Perhaps these will be badly behaved coefficients, so that if 
there is an associated "Hilbert-Polya" Schrodinger operator, its potential will be very badly 
behaved, and cannot be treated function. 

In any case a desire to explain GUE suggests further consideration of the question: 
what conditions must a Schrodinger potential V{x) satisfy in order to have eigenvalues 
that reproduce the Weyl asymptotics of the zeta zeros and also obey local normalized GUE 
statistics? 

6. Concluding remarks 

From the perspective of de Branges spaces, the natural operator to consider is a canonical 
system. Indeed, if the Riemann hypothesis holds, a canonical system "Hilbert-Polya" opera- 
tor must exist with structure function (15. Sh whose (real) eigenvalues are the imaginary parts 
of the zeta zeros. This paper is motivated by the observation that under some conditions a 
canonical system can be (nonlinearly) transformed to the familiar form of a one-dimensional 
Schrodinger operator on a half-line, also known as a singular Sturm-Liouville problem. This 
raises the possibility that the "Hilbert-Polya" operator proposed for zeta zeros might be 
realizable by such a Schrodinger operator, which falls squarely in the framework of quantum 
mechanics. 

The existence of this nonlinear transformation reflects the fact that the de Branges 
theory generalizes part of Weyl's theory of singular Sturm-Liouvlille problems, cf. Everitt 
and Kalf [29]. The Morse potential example considered here is a singular Sturm-Liouville 
problem which falls within Weyl's theory. That theory assigns to each such problem an 
analytic function, the principal Weyl-Titchmarsh m-function, which encodes information 
about the spectral data of the problem, encoded in the initial conditions. This function is 
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discussed in Appendix B. This function is analytic in the upper half-plane and has positive 
imaginary part there, i.e. it is a so-called Herglotz function. 

One may also associate directly to a de Branges structure function E{z) the function 

^E(')--=-J^y (6-1) 

which we will here call the de Branges m-function. It is a Herglotz function associated to 
the de Branges multiplication operator (M^jP^) which plays a role analogous to that of 
the principal Weyl-Titchmarsh m-function for Schrodinger operators. That is, it encodes 
certain information on the spectrum of self-adjoint extensions of this symmetric operator. 
More precisely, under the de Branges transform it equals a "classical" m-function attached 
to a Dirac operator, in the sense of Levitan and Sargsan [45\ Chap. 3]. This Dirac operator 
is directly constructed from the canonical system (j5.2p . by a recipe we omit here. In the 
case of the de Branges structure function ()5.6p above, using a Whittaker function identity 
[H (13.4.30)], we obtain 



Biz) zW.^^{e-) 1 



A{z) VFi.,(e«) 



. (e") 

'( i:!l_ + e--i) 



(6.2) 



For comparison, in Appendix B we determine the principal Weyl-Titchmarsh m-function 
for the Morse potential on the right half-line (Theorem 18. Ih . The formula (j6.2p should be 
compared with Theorem 18.11 with parameter k = The computations in Appendix B 
may be useful as a guide for determining relations between Weyl-Titchmarsh m-functions 
and suitable de Branges m-functions. 

7. Appendix A: Whittaker Functions in the Complex Domain 

We recall here facts about Whittaker functions Wk,ij,{z) as functions of three complex vari- 
ables (k, ^, z). They are a type of confluent hyper geometric function, and are single- valued 
in the variables k and fi, and multi-valued in the z variable, with a singular point at z = 0. 

The standard confluent hyper geometric function iFi{a, (3; z) is given by the power series 
expansion 

T{(3)^T{a + j)zi 



,F,(a,/3,.).= ^^^^,^-. (7.1) 

The normalized function 



iFi(a,/3;z) := ^iFi(a,/3;z) (7.2) 

is an entire function of three complex variables (a,/3, z) € C^, as can be read off from the 
uniform convergence properties of the expansion (j7.ip on compact subsets of C^, observing 
that entire and ^r^^jf^ is a polynomial in a, cf. Buchholz [8, Sec. 1.3]. (The 

unnormalized function iFi(a, /3; z) has simple poles at /3 = 0, —1, —2, ... for most (a, z).) It 
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satisfies the confluent hypergeometric differential equation 

d'^F dF 

which has a regular singular point at z = and an irregular singular point at z = cxd, of 
irregularity index one. (Slater [BJJ). 

The Whittaker functions Mk^^{z) and Wk,,^{z) were given by Whittaker [TSj in 1904 as 
particular solutions to the differential equation 

This equation also has a regular singular point at z = and an irregular singular point at 
z = oo on the Riemann sphere, of irregularity index one. Whittaker sets 

M«,^(z) := e-tz5+/^iFi(^ - k + ^l,2^l + l-z), (7.5) 

with z~^ = e'^'°s^, and the z-plane is cut along the negative real axis. This defines the prin- 
cipal branch of the Whittaker M-function. Under analytic continuation in the z-variable, 
this function is multivalued, with multivaluedness arising entirely from the function z~^. It 
also has poles at = 0, — ^, — 1, ... and to eliminate these (excluding /i = 0) Buchholz [Sj p. 
12] introduces the normalized function A4k,^(z), defined by 

M.,^{z) := ^,^^^_M,,^(z). (7.6) 

The Whittaker function Wk,^(z) is specified by the asymptotic property of having rapid 
decrease as z = x — > oo along the positive real axis. Whittaker [73] defined it using two 
integral representations given in Whittaker and Watson [7i\ Sec. 16.12]. In terms of the 
confluent hypergeometric function above we have (Truesdell \70\ p. 170]) 

W^,^{z) = -^-f^ — iFi (--« + /., 2/x + l;z) 

+ — M.-2. + 1;.). (7.7) 

r(i-K + /i) 2 

This formula can alternatively be written 

W^,,{z) = -^tltL_M^ (^z) + /^^^^ , M,,-,{z), (7.8) 

r(i-K-/u) r(i-K + ^) 

which incidentally shows that WK^^iz) = Wi^^-^{z). The formula (j7.7p directly defines 
Wfc^fiiz) as an analytic function of three complex variables, in the region (k, fi) G and all 
z in the complex plane cut along the negative real axis. ( The analyticity is read off using 
the function (j7.2p .) We term VFk,^(z) on this region the principal branch of the Whittaker 
function. Under analytic continuation of this function in the z-variable, this function is 
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multivalued, and becomes single-valued on a logarithmic covering surface around the point 
z = 0, remaining an entire function of k, and fi on every branch above a fixed z. (For certain 
values of n and /x it may be single-valued on a finite cover of C* = C\{0}.) 

The function W^^^^iz) also has the contour integral representation ([74:, Sec. 16.12]) 

W,,^{z) =r{l + K- /i)e-iz'^ f--^ / + l)K-i+Me-tdA , (7.9) 

2 \ 2711 Jc z J 

in which the contour C goes from +oo to 0^ and back by a "keyhole" contour around the 
positive real axis, keeping t = —z outside the contour, with arg{z) taking the principal 
value, \arg{—t)\ < vr, and arg(l + ^) taking that value which goes to zero at t — > by a 
path outside the contour. If we restrict z to the positive real axis, one obtains by estimating 
the contour integral ()7.9p that 



W^,^{z) = 0[e 

viewed as a function of /i, for fixed k and z, a constant C depending on n and z. This 
implies that Wk,^(-z) is an entire function of order at most 1 in the /i- variable. 

The distinguishing property of the principal branch of W^^^iz), is its rapid decrease 
along the real axis z = x as x ^ oo. All other linearly independent solutions to Whittaker's 
equation ()7.4p increase rapidly (in absolute value) along the positive real axis as x — > oo. The 
rapid decrease property of the principal branch of the Whittaker function holds uniformly 
as \z\ — > oo on the angular sector | arg(z)| < f — e for any fixed e > 0. Furthermore, for 
fixed {K,fi), it has on the larger angular sector — vr + e < arg(z) < vr — e, an asymptotic 
expansion ([TH Sec. 16.3]) which gives 

W^A^) = e'^^z- [\^0 {l^^^l^^ ^ . (7.10) 

The Whittaker function solutions having this rapid decrease property glue together smoothly 
in the (k, /i) variables (uniformly with |Kp -|- < i? for any fixed i?) to form entire func- 
tions of (k, h) G C^, yielding the principal branch defined above. For certain values of k and 
/i we have a finite covering of C*; in these cases there are a finite number of other branches. 



8. Appendix B: Weyl-Titchmarsh m- Functions for Morse Potentials 

For Schrodinger operators on the half-line [uo,oo) with a real- valued potential, the Weyl- 
Titchmarsh m-function m{uQ, E) is an analytic function defined onE ^ C\M, which encodes 
information about the spectrum of the Schrodinger operator for various constant boundary 
conditions. It maps the upper half-plane ^(-B) > into itself, and its values on the lower- 
half plane Im{E) < satisfy m{uo,E) = m{uQ,E). Spectral information is extractible 
from its asymptotics in the upper half plane as the real axis is approached, see Gilbert 
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and Pearson [32j . In general it is not defined on the real axis, but in the case where the 
spectrum is pure discrete then it has a meromorphic extension to the whole plane C, which 
is real- valued on the real axis. 

To define it, let 9a(u;E), (pa{u,E) denote the solutions to the initial value problems 

d2 



+ V{u) ^{u) = E^{u) 



with 



eo.{uo,E) 

0'^{uo,E) 



cos a 
sin a 



and 



fa{uO,E) 



— sma 
cos a 



We define the m-function ma{E) := ma{uQ,E) by the condition that 

ip{uQ,E) := eaiuo,E) +ma{E)ipa{uo,E) 



(8.1) 



belongs to L'^{[uo,oo);du); there is a unique solution when Im{E) ^ 0. The case a = 
gives the principal Weyl-Titchmarsh m-function, where ipo satisfies Dirichlet boundary 
conditions at the left endpoint, and Oq satisfies Neumann boundary conditions; we write 
m{uo,E) := mQ{uQ, E). For the principal m-function we have the alternate formula 

i^'{uo,E) 



m{uo,E) 



d 



ip{uo,E) du 



logilj{u,E) 



U=U() 



.2) 



which follows since ^{uo,E) = 1 and %1)'{uq,E) = m{uQ,E). The different m-functions for 
different a are related by linear fractional transformations, namely 



maiuo,E) 



{cos a)inQ{uQ, E) — sin a 



(8.3) 



{sma)mQ{uQ, E) +cosa 

the principal m-function determines them all. With our definition the m-function maps 
the upper half-plane = {E : '^{E) > 0} into itself, so that it is a Herglotz function. 
(Our definition is the negative of the m-function for a given in Titchmarsh \67\ Sec. 2.1], 
which maps the upper half-plane to the lower half-plane C~, cf. [BT, Sec. 3.1].) As noted 
above, in the general case an m-function need not analytically continue across the real 
axis anywhere. However in the case at hand, that of pure discrete spectrum, the Weyl- 
Titchmarsh m-function extends to a meromorphic function on C, and all of its singularities 
are simple poles, located on the real axis. This case is treated in Chapter II of Titchmarsh 
[57] . The principal m-function m{uo, E) then has a pole at any point where tp{uo, E) would 
be a multiple of 9j(no, E). 

Theorem 8.1. The Schrddinger equation on the half-line [uo,oo) with Morse potential 



Vkiu) 



+ ke^ has principal Weyl-Titchmarsh m-function m{uQ,E) with eigenvalue 



E = for endpoint uq given by 



m(no, z'^ 



"0 1 



w. 



kAz 



(e«o) 



+ 



+ k 



14) 
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Proof. Since the Morse potential on the half-hne [uo,oo) has pure discrete spectrum, the 
m- function wih be a meromorphic function of the variable E = z^. Theorem 12.1( 1) shows 
that for '^{E) > the subordinate solution is given by 

We conclude that 



m{u,z ) 



tp{u, z"^) 



in which W^^i^(x) = ^VF«:,«(a;). We now use the identity ([H 13.4.33]) 

^K,^.{x) = {\^ - ^^)W.A^) - Ty.+i,^(x). (8.6) 

Taking x = e" and substituting this in (jS.Sp results in 

1^, , , l\ W^i-..,^(e") 



the desired formula. 



m{u, z^) = -e" + fc - - - ^""'7 / , (8.7) 



Remarks. (1) The principal m-function contains enough information to uniquely recon- 
struct the potential V{u) on [uq,oo), under very general conditions; this is a particular kind 
of inverse spectral problem. In the case where the Schrodinger operator has pure discrete 
spectrum, this follows from the fact that the Dirichlet problem eigenvalues are the poles of 
m{uQ,E), and the Neumann problem eigenvalues are the zeros of m{uQ,E); necessarily all 
of these fall on the real axis. The Dirichlet and Neumann spectrum together are known to 
be enough to determine the spectrum. 

(2) An important feature of the principal m-function is that it satisfies a Riccati equation 
in the n-variable, 

— m(n, E) + m{u, Ef = V{u) - E. (8.8) 
du 

This can be exploited in inverse spectral reconstruction of the potential. 

(3) It is well-known that the principal m-function attached to a Schrodinger operator 
on a half-line having a potential V{u) that is locally (and which is in the limit point case 
at 00) has strong restrictions on its asymptotics in the upper half-plane as ^(-E) — > 00. Let 
E = z^ with z confined to a sector < e < arg z < ^ — e. Then for a > and < lio < a, 
for sufficiently large \z\ with e < arg(z) < ^ — e, there holds 

f-a-uo / 1 \ 

m{uQ,z^)=iz- j V{t + UQ)e^"^dt + 0[—-A, (8.9) 
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where the constant in the O-symbol depends on a and e. The basic result asserting m{uQ, E) 
equals i^/E plus an error term is due to Atkinson [3l Theorem 2 fF.], who deduced it using 
(jS.Sp . Atkinson used a different definition of m-function, which in terms of the definition 
used here is ^^^-^ ■ (That is, Atkinson used the Neumann boundary condition in defining 
his m-function, while the principal m-function is based on the Dirichlet boundary condition 
at the left endpoint.) A result in the form (|8.9|) appears in Rybkin |58t Corollary 5.2]. 
Combining Theorem 18.11 with the asymptotics (j8.9p provides information on the relative 
sizes of two Whittaker functions Wi-k^^ie^) and H^_fc^^(e") for in the second quadrant 
region of the /i-plane. 
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